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Abstract. Gabriel F. Calvo and Antonio Picon defined a class of operators, for use in 
quantum communication, that allows arbitrary manipulations of the three lowest two- 
dimensional Hermite-Gaussian modes Tir = {|0,0), |1,0), |0, 1)}. Our paper continues 
the study of those operators, and our results fall into two categories. For one, we show 
that the generators of the operators have infinite deficiency indices, and we explicitly 
describe all self-adjoint realizations. And secondly we investigate semiclassical approx- 
imations of the propagators. The basic method is to start from a semiclassical Fourier 
integral operator ansatz and then construct approximate solutions of the corresponding 
evolution equations. In doing so, we give a complete description of the Hamilton flow, 
which in most cases is given by elliptic functions. We find that the semiclassical approx- 
imation behaves well when acting on sufficiently localized initial conditions, for example, 
finite sums of semiclassical Hermite-Gaussian modes, since near the origin the Hamilton 
trajectories trace out the bounded components of elliptic curves. 



I. Introduction 

In a series of papers, Gabriel F. Calvo, Antonio Picon, and co-authors studied the 
manipulation of single-photon states for purposes of quantum communication, from theory 
to experimental design [1], [5], [6]. They demonstrated the limitations of the metaplectic 
operators when acting on spatial transverse-field modes, and they showed how one can 
overcome those limitations with a different family of transformations [5]. In particular, 
they considered the subspace of states spanned by the three lowest Hermite-Gaussian 
(transverse spatial) modes Tir = {|0, 0), |1, 0), |0, 1)} and found a class of operators that 
would allow arbitrary manipulations of such modes. However, their class of operators 
includes "non-Gaussian transformations" , which are not metaplectic operators. Thus the 
question arises of what properties of metaplectic operators may be extended, at least 
partially, to their non-Gaussian transformations. This could lead to the construction of 
the associated optical system [5]. 

The theory of metaplectic operators is often presented in terms of the Stone- von Neu- 
mann theorem (see, for example, the book of G. B. Folland [9]). Rather than attempt to 
extend the Stone-von Neumann theorem to non-Gaussian transformations, in this paper 
we provide an alternative approach, based on the theory of semiclassical Fourier integral 
operators. We find that a Fourier integral operator ansatz provides approximate solutions 
to the evolution equations of Calvo and Picon, and in the process we will determine the 
canonical transformations associated with the non-Gaussian transformations, in the sense 
of Egorov's theorem. In fact, the canonical transformations are given by the Hamilton 
flow of the semiclassical Weyl symbols of Calvo and Picon's [now semiclassical] differen- 
tial operators. The Weyl symbols are dependent on the semiclassical parameter so 
the Hamilton flow is also /i-dependent. We take this point of view essentially for three 
reasons: (1) The operators may be exactly reconstructed from their (/i-dependent) Weyl 
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symbols. (2) The flow leaves invariant a disc of radius ~ v^, which closes up as /?. 0. 
In the context of laser physics, 




where Wq denotes the radius of the laser beam's waist. Outside this disc of radius ~ Wq, 
the flow goes to infinity in finite time. This is because the flow propagates along elliptic 
curves having two components, one bounded and one unbounded. (See Section I VI I and 
the Appendix.) And (3) the appropriate version of Egorov's theorem has an error of order 
0{h'^), rather than the more typical 0{h). (See Section rVIIII ) Moreover, in Section [IXl 
we give examples suggesting that we have creation and propagation of singularities along 
the /i- dependent flow. 

In this paper we find that there are fundamental difficulties with the generators of 
the non-Gaussian transformations. First of all, we show that, although they are clearly 
symmetric, they have infinite deficiency indices. We are however able to describe all 
self-adjoint realizations in terms of certain boundary conditions at infinity. And there 
are difficulties with constructing semiclassical approximations of the unitary groups (the 
non-Gaussian transformations). We can still find approximate solutions to the evolution 
equations, but the associated symplectic transformations can blow up in finite time, as 
they are described in terms of the Weierstrass p-function, which, as is well-known, has 
a double pole in every period. However, we show that the semiclassical propagator be- 
haves reasonably well when acting on sufficiently phase-space localized initial conditions, 
for example, finite sums of semiclassical Her mite- Gaussian modes. This is because, as 
mentioned above, the Hamilton flow leaves invariant a disc of radius ~ Vh centered at 



In Section [TTl we review the non-Gaussian transformations of Calvo and Picon and show 
that the generators have infinite deficiency indices. We then describe all self-adjoint re- 
alizations in terms of certain boundary conditions at infinity. In Section IIIII we introduce 
the semiclassical version of Calvo and Picon's work, so that we can solve their evolution 
equations approximately, in the semiclassical regime. We illustrate this method in Sec- 
tions [IV] and |V] by first taking simpler generators and using semiclassical Fourier integral 
operators to find exact solutions to the associated evolution equations, corresponding to 
"Gaussian transformations". We use the same method for the more difficult non-Gaussian 
transformations in Sections I VI I and IVIIt but in that case we only have approximate solu- 
tions to the evolution equations. In Section IVIIII we discuss Egorov's theorem and give a 
property of non-Gaussian transformations that is analogous to a property of metaplectic 
transformations. We give concluding remarks in Section IIXI and we include an appendix 
for relevant facts regarding elliptic curves. 

Notation: We write = i^r- and similarly for the other variables. Also, we will use 

I ax ' 

the semiclassical Weyl quantization of a symbol a, defined by 



We will mostly be in dimension n = 1 or n = 2, and we will mostly deal with polynomials 
a (hence resulting in semiclassical differential operators with polynomial coefficients). 



the origin. 
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II. The Non- Gaussian Transformations of Calvo and Picon 

In this section we recapitulate the recent work of Calvo and Picon [5]. They introduced 
the following eight generators acting on Hermite-Gaussian modes: 

% = --{aldy - aja^), 



1 

'h = 

1 



— Q (^I ~^ ^-r^x^x Oi^x^x^x ^1^2/^2/ (^x(^y(^y) 1 

i 



% = -^{dl - - d\d\dx + aj,a^a^ - d\dldy + a^.a|^aj^), 

' Q ~^(^CLy I 0/yOjyO/y QjyQjyQ'y CtyCtr^Ct^ CtyCLr^Ct^^ ^ 

^ z 

" ^y ClyCiyCly ~\~ dyClyCly (XyQj^Qix H~ dy^x^^)^ 

% = + 3(a^a^ + dldy)], 

defined in terms of the creation and annihilation operators aj. = ~ ^^nd d^ = 

+ ^) , respectively, and similarly for the y-variable. We note first of all that 71 and 
% are unitarily related by the Fourier transform JF: 

% = oT^o 

and that Tg and % are obtained from and by simply interchanging the variables x 
and y. 

These generators, within the subspace generated by the lowest three Hermite-Gaussian 
modes Hf = {|0, 0), |1, 0), |0, 1)}, obey the SU{3) algebrjfl 

= ifabc'Tc 

(a, 6, c = 1, 2, ... 8), where the only nonvanishing (up to permutations) structure constants 
fabc are given by 

/l23 = Ij /l47 = /l65 = /246 = /257 = /345 = fsiG = 1/2, and /458 = /eTS = 

We note that the triad of generators 

gives a SU{2) group that conserves the mode order. The remaining two SU{2) groups 
are formed by the triads 

and 

Unitary operators IJy-^ generated by the first triad give rise to superpositions between the 
two modes |1, 0) and |0, 1), leaving invariant the fundamental mode |0, 0). Unitarities 11^2 



^We sum over c, which is only relevant here for (a, h) = (4, 5) and (a. h) = (6, 7). 
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and Uts, generated by the second and third triad, produce superpositions between the 
two modes |0, 0) and |1, 0) (leaving invariant |0, 1)), or the modes |0, 0) and |0, 1) (leaving 
invariant |1,0)), respectively. For reference, we state the operations for Ur^' 

,t . . ,t 





0,0) 




0,0) 




0,0) 




1,0) 




1,0) 




1,0) 



'2' 
'2' 



'2' 



|0,0) 



'2' 



e2|l,0). 



'2' 



,t_ 

'2' 



The corresponding operations for U^^ and Uy^ are similar. 

However, a difficulty arises when attempting to extend to larger subspaces of L^(M^): 
starting from the smallest natural domain, the domain of [finite] linear combinations of 
Hermite-Gaussian modes, the symmetric operators 74, Ts, Te, and % all have multiple 
self- adjoint realizations. 

We first prove that the deficiency indices of 7^ are both infinity. For this we consider 
how ^4 acts on the basis of two-dimensional Hermite-Gaussian modes, which we write 
either as |m, n) or as V'm; with m being the mode order in the x- variable and n being the 
mode order in the y-variable. We have that 



(1) 

where 



We take the domain of 74 (as an unbounded linear operator) to be the subspace D{T4^) of 
[finite] linear combinations of Hermite-Gaussian modes. This domain is dense in L^(M^), 
and 71 is clearly symmetric with this domain. And it is easy to check that the domain of 
the adjoint operator is 

{00 00 
9 = E 9X e L\R'y, J2 Wm^^ + /?;;-i^:;-ir < . 
m,n=0 m,n=0 ) 

and that, for g = TZ,n=Ml. e ^(t*), 

00 

m,n=0 

And we occasionally find it convenient to use the formal operator on sequences given by 

£4 : {gTm - iUgZ := PZg:^+, + ^il^.^gl^.^. 
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Two somewhat degenerate cases occur when n 
separately. For n = 0, we get 



and for n = 1, we get 
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and n = 1, so we treat these 
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In all other cases we have 




























V 



■J 



where < for all m and n. 

Hence 71 acts like a Jacobi matrix for any fixed n, fixing the mode order in the y- 
variable, and 7^ may be decomposed accordingly. To be precise, we consider the orthog- 
onal projection operator onto the n*'' mode in y: 



(3) 



m=0 



Then ([T]) shows that 71 is decomposed according to 

where 7i„ and are by definition the kernels of 1 — P„ and respectively. That is, 
we have (see, for example, the book of Kato [T3] ) 

PnD{t)cD{%), tn^cHn, and tU^cH^. 

Now, to prove that is not essentially self-adjoint on D{'T^, for any fixed n we consider 
the equations 



ZUr. 
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for z G C. Tliis is a recursion relation for Mm+i in terms of and Um-i- In the special 
case when n = we simply take Uq = Ui = 0; then the solution is uniquely determined 
by the initial value U2- If n = 1 we simply take 'Uo=0, so that the solution is uniquely 
determined by the initial value ui. In all other cases the solution is completely determined 
by M_i = and the initial value mq. Hence Um is the initial condition multiplied by a 
polynomial in z, which we write as Pm{z). 

Let z G C\M and denote by N- the orthogonal complement of 7^(74 — z/), called the 
deficiency subspace of the operator T^. This subspace is precisely the subspace of solutions 
of the equation 

f*(j) = Z(j). 

Considering D{72), stated in ([2]), we see that is the subspace given by solutions of the 
difference equation 

{Uu)m = ZUm, M_i = 

(or the appropriate modified version for n = and n = 1. For the remainder we always 
write the initial condition as uq, for convenience). Hence for fixed n the deficiency subspace 
is at most one-dimensional; moreover, it is nonzero if and only if 

that is, if 

oo 

\Pm{z)\^ < OO. 

m=0 

Non-self-adjointness is then a consequence of the following theorem about the infinite 
Jacobi matrix 



L = 



[bo 


ao 








...X 


ao 


bi 


ai 










ai 


b2 


a2 










a2 


bs 












■••/ 



where aj > and bj G M for all j. We can reduce to this case if we multiply T4 by — 1. 
This theorem is from Berezanskii's book ([3], p. 507), where one may find a beautiful and 
elementary proof. 



Theorem 1. Assume that \bj\ < C (j=0,l,. . . ), aj^iaj+i < a| beginning with some j , 
and 



00 ^ 

Then the operator L is not self-adjoint 



■ < 00. 

a,- 

j=o ^ 



It is elementary to show that the hypotheses are satisfied for am = —Pm bm = 0, 
for any fixed n. Hence for any fixed n the deficiency indices of the resulting Jacobi matrix 
are both 1, and so, summing over n, the deficiency indices for 71 are both infinity. 
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As for the operator Ts, we may either use a shght modification of the above methods, 
or we may simply use the fact that the Fourier transform intertwines 7^ and 7^. Then, 
simply by interchanging the roles of the variables x and y, we see that the operators 
and % also have infinite deficiency indices. 

The symmetric cubic operators ^4, Ts, and of Calvo and Picon do however admit 
self-adjoint extensions. Moreover, using the results of Allahverdiev [2j we may explicitly 
classify all of the self-adjoint extensions in terms of certain boundary conditions at infinity. 
We begin by studying the restriction of 74 to the subspace Tin given by the projection P„ 
as in ([3]). We simplify notation by then omitting "n"; in particular, ipm now denotes the 
(m, n)*^ Hermite-Gaussian mode, and 13m will also lose the superscript "n" . 

For 

00 

m=0 

and 

00 

we denote by [g, h] the sequence with components [g,h]m given by 
We then have Green's formula: 

M 

'^{ihg)mhm - gmihh)m} = "[fi-, 

Since the sequences {g)m-, {h)m, {iAg)m, and {t4h)m are all in £^(No), we then have that 
the limit 

[9,h]oo = lim [g,h]M 

exists and is finite. Hence 

iUglh) - {g\Uh) = -[g,h]^. 
Now (still for a fixed n) we let 

u = [um) and v = (Vm) 

be the solutions of 

Pm-iym-l + Pmym+1 = (m > 1) 

satisfying the boundary conditions 

Mo = 1, = 0, fo = 0, and Vi = —. 

(We make the appropriate trivial modifications for n = and n = 1.) We have that 
u,v e D{f^) n Hn] in fact, 

ihu)m = for all m, (t4t')o = 1, and {iiv)m = 0, m > 1. 
With this set-up we have the following results of Allahverdiev |2] : 
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Theorem 2. (Allahverdiev [2J.) The domain of the closure of T4 restricted to Ti. 
consists precisely of those f G D{T^) fl 7i„ satisfying the boundary conditions 

[/,m]oo = [/,f]oo = 0. 



•n 



For f e D{T^)nnnwe now define 



ri/ 



00 



and 



^2f = [f,u]oo. 

Then, in the precise sense of unbounded operators on a Hilbert space, 

Theorem 3. (Allahverdiev [2J.) The operatorsTi, T2 are (complex-valued, symmetric, 



And now that we have appropriate boundary values, we have the following well-known 
description of all self-adjoint extensions: 

Theorem 4. (Allahverdiev [2].) When restricted to Tin, every self-adjoint extension 



for G M U {00}. For h = 00 the condition ^) should be replaced by [/, m]oo = 0. 
Conversely, for an arbitrary h E M. U {00}, the boundary condition ^ determines a 
self-adjoint extension on Tin- 

Remark 1. Allahverdiev [2] additionally considers maximal dissipative and accretive ex- 
tensions of infinite Jacobi matrices. These correspond to h E C such that Imh > and 
Imh < 0, respectively. He also gives applications to scattering theory. 

It remains to consider 12 as acting on the entire space 



linearly independent) boundary values of 71 restricted to 7i. 




00 



n=0 



Let denote the restriction of 7^ to 7-^„, that is, the operator in Tin with D{T^) = 
DiT^ n Tin such that T^f = 7lf G 7i„. As we have just shown, the closure of is 
obtained by extending to the space of / G D{T^*) satisfying F"/ = V^f = 0. We now 
prove the analogous result for the larger space. 
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Theorem 5. The closure ofT^ is obtained by extending to the space 
Vci = {fe D{f:)- Vl{PJ) = mPJ) = Wn}. 

Proof. The closure of T4 is of course the adjoint of 7^*, so we are to show that V^i is equal 
to 

D{f;*) = {Fe L\R^y, 3H G L\R^) such that {f;g\F) = {g\H)^g e D{f;)}. 

Moreover, we note that D{Tl*) C D{T^) since T^* and 7^ have the same adjoint and 
since T^* is automatically symmetric. 
Now let F e72 and write 

F = Y,K^'L 9 = Y.9':nrn., and// = 5^ 

Then we wish to find those precise conditions on F such that there is some H G L^(M^) 
with the property that 

{f:g\F) = {g\H) Wg G D(t;). 
Simply by restricting to g E DiT^j, we see that it is necessary to have 

TTTi _ on -pn _|_ on -pn 

So we must find the conditions on F such that 

{f:g\F) = {g\f:F) Vg G D(t;). 

We recall that 

M 



SO we see that 



JM 
n=0 

converges absolutely. Moreover, by the dominated convergence theorem. 



n=0 n=0 n=0 

Hence 

00 

Y,[9,Fr^ = (glVF) - {f:g\F) Wg G D(t;). 

n=0 

So D{f^*) is precisely the set of F G D{f^) such that 

OD 

E[^'^]- = o v^7GD(t;). 

n=0 

But then this is equivalent to 

[g, F]^ ^ [^7, Pn(F)]^ = Vn, V^? G /^(t/) n H., 
which in turn, as shown by AUahverdiev [2], is equivalent to 

r^(P„F) = r^(P„F) = Vn. 
So indeed we have V^i = D{f;*). □ 
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We have also seen that self-adjoint extensions of T^^ are in one-to-one correspondence 
with boundary conditions of the form 

r^(/) - KTlU) = 

for /i„ G M U {c>o}. Now we define 

= {fe D{f:)- ri{pj) - /i„r^(p„/) = vn}, 

and we denote the extension of to this domain as T^. We next show that all self-adjoint 
extensions of 7^ are of this form. 

Theorem 6. Every self-adjoint extension ofTl is determined by extending the domain 
to a set of the form 

Vh = {fe D{f:)- TUPnf) - h^T^iPJ) = Wn}, 
where h = {hn)^=Q is an arbitrary sequence with /i„ G M U {c>o}, and by the rule 

= {uf) ^ f:f. 

For hn = oo the condition should be replaced by (Pn/) = 0. Conversely, for an arbitrary 
sequence h = {hn)'^=Q with /i„ G M U {oo}, the set Vh determines a self-adjoint extension 
of%. 

Proof. We begin by proving the converse: we take an arbitrary sequence h with hn G 
M U {oo} and prove that the extension to the domain T>h is self-adjoint. 
We first show that Vh is symmetric, that is, we show that 

{f:g\f) = {g\f:f) yf,geVh. 

As in the proof of Theorem 0, we use the identity 

oo 

J2i9jroo = {9\vf)-{f:g\f). 

n=0 

But now we may use the identity 

[9, nil = [g, «1m[7, vTm - [g. vTmU. uTm. 

where and are the functions occurring in the definitions of F" and F2. Hence, in 
the limit, 

[gJToo = [TiPngW^Pj) - iv^PngmPuD- 

Since / are in D^, we can use the identities 

niPnf) - h^TliPJ) = 

and 

n{Png) - hj:i{Png) = 
to see that = for all n. Hence is a symmetric operator. 

To show that is a closed operator, we suppose that 

VhBfk^f in 

and 

f^fk^F inL2(M2). 
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Then clearly, for all n, 

Pnfk Pnf in Hn 

and 

PnXtfk = Xt^rifk ^ PnF in Hn- 

But since 7^^'", the restriction of Tj^ to Tin, is closed (as is T'^), we then have that 
Pnf eVhOHn (hence / G Vh) and that fi'Pnf = Pnf if = PnF (hence %^f = F). So 
7^'^ is a closed operator. 

Next we show that the deficiency indices of Xt are both zero. Suppose u± G D{T^*) are 
such that u± G Ker('7^^* =Fi). Then PnM± G Ker('7^'*'"'* ^i), and, since T^'^'" is self-adjoint, 
PnU± = for all n. Hence u± = 0. So the deficiency indices of Tj^ are both zero, and Xt 
is a closed operator. Hence Xt is self-adjoint. 

For the other direction, let T^, with domain T>e, be some self-adjoint extension of T^, 



-A 

v,nnn3 fk^ f in 7^„ 



and let T^*^'*^ be its restriction to Hn- To show that T^'"' is closed, we suppose that 



and that 

Then, in particular, 
and 



T^' fk-^ F in Hn- 
V,3fk-^f in L^(M2) 
f:fk ^F in L2(M2). 



Since is closed, we have f E H Hn and '7^^'"/ = F, which proves that T^'"' is a 
closed operator. And since the deficiency indices of Xt are both zero, it is clear that the 
deficiency indices of X^'" are both zero. Hence X^'"' is self-adjoint. 

Now 7^^'" is a self-adjoint extension of X^, so the results of AUahverdiev cited above 
show that X^'^ must be given by a boundary condition of the form 



r" (/) - hnV^M) = 

for some G M U {oo}. Hence we have C T>h, that is, 

j-e ^ q-h 

J-A -'4 ; 



and since both operators are self-adjoint, we then have X^ = X^. This completes the 
proof of the theorem. □ 



We have now categorized all self-adjoint extensions of the operator 74, so we are pre- 
sented with the basic question: which is the "right" extension? We expect the physics of 
the problem to dictate the appropriate extension, for which we may return to the origi- 
nal work of Calvo and Picon [1],[5],[6]. They state that the cubic generators, 7^, 7^, 7^, 
and 77, can be implemented with passive optical elements having higher-than-first-order 
aberrations (nonquadratic refractive surfaces) [5]. Perhaps the physics of the apparatus 
will determine the "right" self-adjoint extension. 
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III. The Semiclassical Formalism 

In the case of Gaussian transformations, one may use the theory of metaplectic oper- 
ators (sometimes presented in terms of the Stone-von Neumann theorem) to deduce the 
underlying canonical transformations. However, non-Gaussian transformations are not 
metaplectic operators, so the arguments must be modified. Calvo and Picon then ask if 
the Stone-von Neumann theorem can be extended to the case of the above cubic genera- 
tors, for then "this would enable one to find the explicit form of the symplectic transform 
and thus the construction of the associated optical system" [5]. In the following sec- 
tions we take a different approach: we approximate the non-Gaussian transformations by 
semiclassical Fourier integral operators, by starting with the semiclassical Fourier integral 
operator ansatz and then by solving the resulting eikonal equation and transport equa- 
tions. This, along with Egorov's theorem (stated in Section IVlIip . justifies the claim that 
the underlying canonical transformation is the Hamilton flow associated to the evolution 
equation. As a first step, in this section we put the work of Calvo and Picon in the 
framework of semiclassical analysis. 

We wish to construct approximate solutions to the evolution equations associated to 
the cubic generators 74, Ts, Tg, and approximate in the sense that we will work in 
the semiclassical regime. For this we start with the two-dimensional semiclassical ground 
state 

|0,0) = (7r/i)-^e^2k(-'+2^'). 
To get the other semiclassical Hermite functions, we apply the creation operators 




We also have the corresponding annihilation operators 




so that 

[a^,aj,] = 1 and [fly, a])] = 1- 
The [normalized] semiclassical Hermite functions are then given by 

|m,n) = (m!n!)-5ai"aj"|0,0). 

We again let 

2^ (^ic CLj.CL^CLx ^x^^^^ dj^dyOiy CixdyCly^ 

but now where the creation and annihilation operators are semiclassical. In this context 
it is more convenient to introduce the semiclassical differential operator 

For reference, we note that 

(6) = \x {{hD^f + [hDyf + x^ + y^) - hx - hh'D,, 
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and we note that the semiclassical Weyl symbol of T4 is 

1 5 
P^ix, y, ^, r?; h) = -x{x^ + y'^ + ^'^ + rf) - -hx. 

We then have the permutation properties similar to those in Section [III Explicitly, for 
"we have 





|0,0) 


= cos (2-V2/,i/2t)|o,o) 


— i sin (2 


-V2/,i/2^)|1^0) 




|0,0) 


= cos (2"i/2/,i/2^)|o,o) 


- sin (2~ 


1/2/, 1/2^)1 1,0) 




|0,0) 










|i,o) 


= cos (2-V2/,i/2i) |i^o) 


— i sin (2 


-1/2/^1/2^)10,0) 




|i,o) 


= cos {2-'H'/h)\l,0) 


+ sin (2" 


1/2/^1/2^)10,0) 




|i,o) 









The operations for f/n and [/pa are similar. Hence, in this choice of scale, we have 
permutations of the three lowest Hermite- Gaussian modes when 

TT 

V2h 



IV. Warm-Up: The FIO Representation of a Gaussian Transformation 

For non-Gaussian transformations, the approximate representation by Fourier integral 
operators will be somewhat complicated, so we begin with a simpler situation: the case 
of a Gaussian transformation. For the sake of concreteness, we restrict our attention to 
the semiclassical differential operator 



92 

P(h) = -h{a\al + d^ciy) = -^^-q^ - ^V^ 

though all ten generators of U{S), the metaplectic representation of S'p(4, M) (see [5j^ 
may be treated in the same way. Our goal is then to solve the evolution equation 



{hDt + P)u = 
u\t=o = V. 



(7) 

Following the method outlined in the book of Grigis and Sjostrand ([10], p. 129), we try 
u = Utv{x,y) = (27r/i)-2 f f e^^(*'^'^'«''')a(t){;(^, r/) d^t^r]. 



where 99 is a quadratic form in (x, y, ^, rj). Here v denotes the semiclassical Fourier trans- 
form: 

'^^{^,V)= / e~T^'^'-^^'^y'^^v{x,y) dxdy. 
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With this ansatz, we arrive at the expression 



{hDt + h^D.^Dy - xy)u = (27ihy 



dip dip dip 
dt dx dy 



-ih{2TTh)~^ jj 

Thus we wish to solve both the eikonal equation, 
(8) 

and also the transport equation, 
(9) 



da d'^ip 
dt dxdy 



eh a{t)v{^,7]) dr] 
e^v{^, rj) d^ drj. 



dip dip dip Q 
dt dx dy 



da d'^ip 

+ ^r-^a = 0. 



dt dxdy 

It is due to the special form of the operator P that we are able to solve this in such a way 
that the amplitude a depends only on t. Moreover, we want ip to satisfy 



and we want a to solve 



f{0,x,y,^,r]) = x^ + yr], 



a(0) 



so that the initial condition is satisfied in the evolution equation ([7]). 

We begin with a solution of the eikonal equation ([8]). The semiclassical Weyl symbol 
of P is 

p{x,y,tv) =^V-xy, 
and we write the semiclassical Weyl symbol of the evolution equation ([7]) as 

9{t, X, y, T,^,ri) = T + p{x, y, ^, r]), 

which is actually independent of t. The eikonal equation may then be written in the 
simple form 



(10) 



g(t, X, y, dtip, d^ip, dyip) = 0, 



whose solution, which we now sketch, is provided by Hamilton- Jacobi theory. 
Hamilton's equations for g are 
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Hence it is natural to identify the evolution parameter with time t. These equations may 
be easily solved to give the Hamilton flow of g: 

x{t) — xo cosh i + 770 sinh t 
y{t) — Ho cosh i + ^0 sinh t 
' r(t) = To 
C i't) — Vo sinh t + ^0 cosh t 
r]{t) — xq sinh t + rjo cosh t. 

One may think of this as just being the Hamilton flow of p, since the flow in the (t, r) 
variables is trivial. In this point of view, we may write the Hamilton flow of p in the 
matrix formulation: 



/xit)\ 




/ cosh t 








sinht\ 


f Xo\ 


yif) 







coshi 


sinht 





yo 


m 







sinht 


cosht 











ysinht 








cosh t J 


\VoJ 



However, for the time being we take the point of view of the Hamilton flow of g. This 
Hamiltonian system (in a six- dimensional cotangent space) is completely integrable, since 
we have the three (= |) independent conserved quantities 



The fact that we have three conserved quantities corresponds to the fact that the flow is 
constrained to a Lagrangian submanifold; that is, a three-dimensional submanifold A of 
the cotangent space such that the restriction of the symplectic form cr to A is zero (i.e., 
is isotropic): 

(t\a = 0. 

To construct a solution of the eikonal equation, the basic idea is to start with a two- 
dimensional isotropic submanifold A' and then to propagate A' with respect to the Hamil- 
ton flow of g, filling out the whole Lagrangian submanifold A. To this end, we let 

A' = {(0, Xo, I/O, To, Co, Vo); 9(0, Xo, yo, To, Co, ^0) = To + CoVo - XoVo = 0} 

= {(0, Xo, yo, xoyo - Co^o, Co, Vo)}, 

which we think of as the "submanifold of initial conditions" , and where we consider (Co, ''70) 
as universally fixed. Propagating along the Hamilton flow of g, we thus get the whole 
manifold A: 

A = {{t,x{t),y{t),x{t)y{t) ^{t)7]{t),^{t),7]{t))}. 

On the one hand, we may think of a trajectory along the Hamilton flow as being 
determined by the parameters (t,Xo,yo)- On the other hand, we may think of it is as 
determined by the parameters {t,x,y), since 

Xo = X scch t — rjo tanh t 
yo = y sech t — Co tanh t. 
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Hence, rewriting the variables {C,o,tio) as {C,,i]), we may rewrite A as 

A = {(t, X, y, (xy — ^r]) sech — (x^ + yq) tanht sech t, 

ytanht + ^secht, xtanht + r^secht)}. 

To conclude the solution of the eikonal equation, we seek a function ip such that A is 
the graph of the gradient of ip, to agree with ( ITOl) . and such that ip{0, x, y, ^, t]) = xC,+ yrj. 
This is easily accomplished, and we thus have the phase: 

ip(t, x,y,^,ri) = {xy — ^r]) tanh t + (x^ + yrj) sech t. 

And one may now check directly that this is the solution of the eikonal equation. 
It is now easy to solve the transport equation ([9]): 

a{t) = secht. 

So we have the following expression for the solution operator Uf. 

Utv{x,y) = secht / / exp I — [{xy — C,!]) tanht + {x^ + yrj) sech t] ) ■0(^,77) 



We may use the method of stationary phase (which, in this case, is exact) to simplify this 
expression and get an integral in terms of v. We withhold the details for now, since this 
will be accomplished in the next section for a different but similar operator. 

Also, it is known from the general theory (see, for example, [10]) that is a generating 
function of the canonical transformation, which in this case is the Hamilton flow of the 
symbol p. That is, the Hamilton flow of p is given by 

dip du) ^ \ ( dip dip 
ot, or] J \ ox ay 

This can also be checked directly, now that we have an explicit expression for the phase 
ip. 



V. The Gyrator Transform 
The exact same method may be applied to the semiclassical differential operator 

Ti = h(dldy + d^yd^) = xy — h? 



dxdy 

In the previous section, the solution of the evolution equation was exact, so the semiclas- 
sical parameter h ultimately played no role. Hence in this section we simply let /i = 1. In 
this case, the solution to the evolution equation 

= 




for t G (— 7r/2, 7r/2) is given by 

(11) u{t,x,y) = {27i)~'^\sect\ jj exp (i[{x^ + yr]) sect — (xy + ^r]) tan t]) v(^,r]) d^ di]. 

It remains to extend the solution to all t G M. 

We may use the method of stationary phase (which is exact in this case) to rewrite this 
integral as a double integral involving v. This simply amounts to a use of the following 
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fact. Let Q be a real, non-degenerate nxn symmetric matrix. Then the Fourier transform 
acts as follows (for details, see [ID], p. 21): 

jr. ef(-'O-) ^ (27r)te'f^s'^Q|detQ|-^e-^<«'«"'«>. 
After some calculation, for t G (0,7r/2) we arrive at the integral expression 

(12) «(t,x,y) = (2vr)-|sint|- J j vja, h) exp + ' + ) rfarffe. 

The right-hand side is known as "the gyrator transform" of v [15]. The benefit of this 
expression is that we may now extend the solution u{t, x, y) to t G (0, n). For t G (7r/2, tt) 
we may again use the method of stationary phase to return to the expression ffTTj) . There 
is no phase shift, since in this example sgnQ = 0. 

Hence we have completely determined u{t,x,y): for t G M\{(2A; + l)vr/2; A; G Z} it is 
given by flTT]) . and for t G ]R\{A;7r; A; G Z} it is given by f[T^ . This is analogous to the 
standard parametrization of the circle by four charts of graph coordinates. In fact, it is 
not only analogous, but intimately related; the exceptional points for (fTT!) (resp. (fT2|) ) are 
precisely those for which the Lagrangian submanifolds, swept out by the Hamilton flow 
of Ti, have degenerate projections onto the {^,r]) plane (resp. {x,y) plane). (For more on 
this phenomenon, one may consult Duistermaat's beautiful review article [7].) 

This unitary group has an important application when t = ±7r/4: it takes Hermite- 
Gaussian modes to Laguerre-Gaussian modes. For this we recall the definitions of the 
extended Wigner transform 

and of the (renormalized) partial Fourier transform 



^2F{x,y) = 




F{x,p)dp. 



Then one may check that 

u{TT/A,x,y) = W{T2v){x, -y). 

In the special case when the initial condition is ^; = hmn, the (m, n)*^ Hermite-Gaussian 
mode, we have 

«(7r/4, X, y) = {-iYW{Km){x,y). 
Moreover, by taking complex conjugates, we have 

u{-Tx/A,x,y) = i''W{h^n){x,y). 

And we recall that W{hmn) is precisely the {m,nY^ Laguerre-Gaussian mode [TB] . 

VI. The FIO Representation of a Non- Gaussian Transformation 

We now turn to the more difficult non-Gaussian transformations. The four non- 
Gaussian transformations used by Calvo and Picon may all be treated by the same meth- 
ods, so we focus on the operator 
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where the creation and annihilation operators are semiclassical. Moreover, we will make 
a slight simplification in order to remove inessential complications. Since the operator 
acts very simply in the y variable, we may instead study the operator 

7q — (fl^ -|- Qjx QjIj^Q^^CIx Oj\flxOjx) 

acting on functions of x only. The following arguments may be repeated for 7^, but with 
some slight changes as outlined in Section IVIII 

We renormalize Tq in order to have a semiclassical differential operator: 

1 3 1 

= -x{x^ + (hD^)'^) - -hx - -ih'^D^. 

The method we use to treat this operator is the same as in the previous sections, but 
there are some complications. The difficulties may be treated by the general theory, but 
the solution is not as explicit and exact. Here we wish to remain in the semiclassical set- 
ting, and we only expect an asymptotic solution to the corresponding evolution equation: 



(13) 



u\t=o = V. 



We will allow the initial conditions v to depend on h. Moreover, for bounded times 
Duhamel's Principle shows that the semiclassical propagator differs from the exact unitary 
propagator by 0{h°°) (after choosing a self-adjoint realization of the generator). 
The semiclassical Weyl symbol of Tq is 

1 3 
Po{x, ^; h) = -x{x^ + ^^) - -hx. 



so then Hamilton's equations are 
(14) 



X = xC, 



Here we have the conserved quantity 

Co =poix,^;h). 

Suppose first that Co ^ 0, so that, in particular, x{t) ^ for all t. Letting 

1 

w = -, 

X 

we have from f|T^ the following differential equation for w: 
(15) {wf = 2CoW^ + 3hw^ - 1. 

If Co 7^ 0, this is essentially the same as the differential equation for the famous Weierstrass 
p-function: 

[p\z)]' = 4[p{z)]'-gM^)-93- 
The two constants g2 and are the so-called "invariants" . We may then solve Hamilton's 
equations, giving the Hamilton flow in terms of the Weierstrass p-function. 
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Figure 1. Hamilton flow lines in the (x,^) upper half plane in the case 
h = 1/10, with variable Co- The lower half plane is obtained by symmetry. 



To be explicit, for Co 7^ we have 
(16) x{t) 



where to is either an arbitrary real constant or an arbitrary real constant plus ^ui, the 
purely imaginary half-period of p (see Appendix). Here p is the Weierstrass p-function 
associated to the invariants 

g2 = ^h^ and g3 



_ l/,3 
4^0 • 



We then also have 



(17) 



x(t) 

pit + to) - \h' 



When /i = and Cq 7^ 0, ^(t) is always strictly decreasing, which follows simply from 
Hamilton's equations (fT4|) . However, when /i > we have more complicated behavior, as 

I Wo, where Wo is in practice the 



shown in Figure [H There is a pocket of radius v3/i 
radius of the laser beam's waist. 



We have very different behavior when Cq = 0. Depending on the initial conditions, we 
have one of the three following cases: 

'x{t) = Vt 

iit) = v^coth (^^{t + to] 
'x{t)=0 Vt 

^(t) = Vshtanh (^^{t + to] 
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or 

'x(t) = iv^sech {V3h{t + to)) 
C{t) = -VShtanh {V^it + toj 

And we have the stationary points (x,^) = (0, ±-\/37i) and (x, ^) = (±\//i, 0). 



Now we look for a solution of the equation f|T3|) of the form 
(18) u = Utv{x) = {2'Kh)-^ j e^^(*'^'«)a(t, x, h)v{^) 

Here we are using the semiclassical Fourier transform, given by 

■O(^) = / e~T^''^v{x) dx. 



We will even allow the phase (p to depend on /i, since we prefer to study the Weyl symbol 
of 2^0) rather than the principal symbol only. 

We take to be a solution of the now /i-dependent eikonal equation 

(10) 'iM^'i,-)-"^- 

We may treat this using the same method as before, using Hamilton- Jacobi theory (see, 
for example. Theorem 5.5 of [10]): given any (xo,^o) ^ I^^, there exists a real-valued 
smooth function ip{t, x; h) (with (^o considered as a parameter) defined in a neighborhood 
of (0, xo) such that (fT9l) is solved and such that 



ip{0,x) = x^o, 
-( 



dip 1 / 2 3 

— (0,xo) = --xo(xo + 4o) + g^-^O' 



^^-(0,xo)=eo. 

However we have not yet managed to get a closed-form expression for the Hamilton flow 
(solely in terms of the initial conditions xq and ^o) or for the phase p. Perhaps this is 
best suited for Laurent expansion methods and numerical methods, but we leave open 
the possibility that one may find closed-form solutions by brute force. 

We next construct the amplitude a. For ffTSl) to be a solution of the equation f|T3l) . the 
amplitude a must satisfy 

ihDt + %){e'^a)=Oih°°), 

which is equivalent to 



And since <y9 is a solution of the eikonal equation, we then have 

dip 
dx 



(20) ( hDt + e-^%oe^ -Po[x,^;h] ]a = 0{h° 
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We simply compute 

e " Toe " -po a;, — ; h 
\ ox 

If we let 

V ■=x— — 

dx dx 

(which we recall depends on /i), then ( l20l) becomes 
It is then natural to define 

at 2 

and 

2i dx dx 
First considering a finite sum, we have 

TV N 

{Ri + hR2) ^ K'an = Riao + ^ h" [Rian + i^aOn-i] + /i^+^i?2a7v, 

n=0 n=l 



with the initial conditions 



so we wish to solve the transport equations 

RittQ = 

Rlttn + R2Cln-l = 0, n>\ 

ao(0,x,O = 1 
an(0,x,O = 0, n > 1. 

There is an elegant solution to the equation 
(21) (^^ + ^divVi + ri^ao = 

which gives interesting geometric information (see j8j and Appendix A of [H]). It is due 
to the fact that 

V = ^(x —-h] — 
\ ' (9x ' ) dx' 



We let 



and we denote the (/i-dependent) Hamiltonian fiow of by k (generated by the Hamil- 
tonian vector field H^^), considered as taking 
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We then note that 



d 

ds 



s=0 



f = vj 



for / e C°°. 
Considering 

ao(t, x,C,)\dx\^^'^ 
as a half-density on A^^^, fl2T|) then becomes 



9t 



+ Cv^ {ao\dx\^^^) = 



where Cvt denotes the Lie derivative. That is, the amphtude, interpreted as a half-density, 
is invariant under the flow. This is the same as 



<(ao(t,x,OMx|i/2 



\dx 



1/2 



Hence 



where dKt denotes the differential of the mapping /tj. 

Now we recall that the mapping Kt may be described in terms of its generating function 
SO that its inverse is given by 

'dip 



Hence 



So finally we see that 



ay 
dxd^ 



{t,x,i\h),i . 



\ dxdi J 



ao(t,a;,0 



Remark 2. The same argument works for more general operators [8], [H]. In particular, 
the simple operators in Sections IIVI and |V] fit into this framework. In Section IIVI we had 
the phase 

Vit, x,y,^,v) = i^y - ^V) tanh t + {x^ + yrj) sech t, 

so that 

det 12 ^ = sech\ 

\ dyd( dydr] / 

And in Section |V] we had the phase 

Vit, X, y, ^, v) = i^^ + yv) sec t- {xy + ^r]) tan t, 

so that 

det 15 ^ Usec^t. 

\ dyd^ dydr] / 
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We may also solve the higher-order transport equations to construct the full amplitude 

oo 

a ~ /i"a„ 

n=0 

in the sense of Borel summation. 

To rigorously solve (fT3|) one must control the error. For this one uses cutoff functions 
in phase space and restricts to initial conditions v that are appropriately localized in 
phase space (for a textbook presentation, see ^). For example, finite sums of Hermite- 
Gaussian modes are localized to the origin in phase space, and so are localized to the 
flow-invariant disc. So fortunately we have good behavior for the physically most relevant 
initial conditions. 



VII. The Operator T4 



For reference, in this section we give the Hamilton flow for the operator T4. The 
methods of the previous section may be apphed, with some slight modifications due to 
the additional parameters. 

The semiclassical Weyl symbol of T4 is 



P'i{x,y,^,r]; h) 
so then Hamilton's equations are 

X = xS, 

(22) 



1 5 
-x{x^ + 1/^ + + rf ) - -hx, 



y = XT] 



KV = -xy- 

Here we have the two conserved quantities 



Co =P4{x,y,^,v) 

ci 



y^ + rf-. 



(23) 



For Co 7^ we have the solution 

1 



x(t) = -Co(p(t + to) + Y^(Cf-5/i) 



where to is either an arbitrary real constant or an arbitrary real constant plus the 
purely imaginary half-period of p (see Appendix). Here p is the Weierstrass function 
associated to the invariants 

^-XCl-hhf and g^ = ^{Cl-hKf ^\cl. 



92 

We then also have 
(24) 



12 



216 



-Pit + to) 



pit + to) + ^iCf-5h) 
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As for y, we have the equation 

1 d ri \ _ ^ 

where x is as in fl23|) . The solution is given by 



and hence 



y{t) = yocos / x{s) ds + rjo sin / x{s) ds 
Jo Jo 

?7(t) = — I/O sin / x{s) ds + rjQ cos / x{s) ds 
Jo Jo 



where of course 

.2 I „2 ^2 



z I z /^z 

We have very different behavior when Co = 0, which breaks into multiple cases. If 
Cf > 5h, then x(t) = for all t. We then have y{t) = and rj{t) = rjQ. As for ^, if 
> 5h, then 

m = -^Cf-5h tan Q^C2-5/i (t + to)) , 

where to is an arbitrary constant. And if Cf = 5h, then either ^ = or ^(t) = ^q^. 

On the other hand, if Cf < 5h, we have one of the following three cases, depending on 
the initial conditions. Either 



x{t) = ±^5h-Cf sech (^^5h-Cf (t + to) 

e(t) = -^5h-Cf tanh (^yju^^l (t + to)) 

y{t) = yocos / x{s)ds + riosm / x{s) ds 
Jo Jo 

?7(t) = — I/O sin / x(s) (is + ?7o cos / x{s)ds, 
Jo Jo 

or, in the case where xo = 0, 

xit) = 0, y{t) = yo, r]{t) = r/o, and 

e(t) = ^Jbh-Cf tanh Q^5/i-C2 (t + to)) , 

or 

x{t) = 0, y{t) = yo, r]{t) = r]o, and 
^(t) = ^hh-Cl coth Q^5/i-C2 (t + t, 



And we have the stationary points 



(a;,|/,e,r?) = (±A/^/i,0,0,0) 



and, for y^ + % = C( < 5h, 



{x,y,^,r]) = (0,?/o,±^5/i - Cf,r]o). 
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VIII. Egorov's Theorem 

When dealing with the semiclassical quantization of a quadratic polynomial, one may 
use the metaplectic representation. This is the method discussed, for example, in the 
work of Calvo and Picon [S]. Here we will simply define the metaplectic representation 
and note the connection with Egorov's theorem. 

Following Folland [9J, we write the Schrodinger representation of the Heisenberg group 

as 

and we write the metaplectic representation as /x. The metaplectic representation is 
sometimes described in the language of representation theory as follows. Let T denote 
the group of automorphisms of the Heisenberg group H„ that leave the center pointwise 
fixed. If T e T, then the composition p o T is a new irreducible unitary representation 
of the Heisenberg group, nontrivial on the center, so the Stone- von Neumann theorem [S] 
says that there exists a unitary operator /i(T) on L^(M"') such that 

^l{T)p{x)^{T)-^ = p o T(x), X G 

Let us compare this to (one version of) Egorov's Theorem: 

e^'*^/'^Opr(g)e-^**^/^ = Oplflg,) 

with qi = q o Kt + 0{h'^). Since we are dealing with unitary operators and the Weyl 
quantization, we have an error of order 0{h'^) rather than the more usual 0{h) (see 
Appendix A of [llj or Section 2 of This is in fact one of the main reasons for our 

using /i-dependent canonical transformations. Here = exp(t_fZp^) is the (/i- dependent) 
Hamilton flow associated to T4. Egorov's theorem is a way of justifying the intuition 
that "the Fourier integral operator e**^/'* quantizes the Hamilton flow of if" . And there 
is a wider class of Fourier integral operators that may be considered as quantizations of 
canonical transformations |lUj . 

This has a useful consequence for optics, where the Wigner transform is a widely used 
tool for studying phase space properties of functions. For this we define the standard 
n-dimensional semiclassical Wigner transform by 



W{f, g){x, = y e f (^x + ^pj g(^x- ^pj dp. 

One may check that for any symbol a one has 

{OpY{cr)f\9) = llcT{x,OW{f,g){x,0 dxdi. 

Then, for example, using Ut, the semiclassical approximate propagator for T4 (with Weyl 
symbol p4, having the Hamilton flow Kj), one has 

{OpYi^WtflUtg) = {e^Op]^ {a)e-^f\g) + 0{h^) 
= {0p^{^t)f\9) + O{h-), 
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where at = cr o -\- 0{h'^). Hence 

jj a{x,OW{Utf,Utg){x,Odxd^ = jj a{Kt{x,0)WU.g){^.i) dxdi + 0{h^) 

= jj a{x,0WU,g){^^^t{x,0)dxdi + O{h^). 
Since this is true for all symbols a, we thus have 

WiUtf, Utg) = wu. g) o ^-t + o{h^). 

If we were dealing with a metaplectic operator, the transformation k would then be a 
linear symplectic transformation, and the identity would be exact. 

IX. Additional Remarks 

Remark 3. So far we have concerned ourselves with the manipulation of Hermite- 
Gaussian modes: those two-dimensional modes generated by applying the x and y creation 
operators d\. and to the fundamental mode 7r~^/^e~2(^ +y ) (one may alternatively use 
the semi classical version). But one may also just as easily consider the manipulation of 
Laguerre- Gaussian modes, generated by applying the creation operators 

i^ = 2-V2(at±,at) 

to the fundamental mode. It is no more difficult to use the above methods, since the 
extended Wigner transform intertwines the two classes of creation operators [16] . 



Remark 4. We have seen that the Hamilton flows of T4, X5, Tg, and T7 blow up in finite 
time. One may now ask: what are the consequences for the propagators themselves? As 
noted by Maciej Zworski [18], the propagators (the non-Gaussian transformations) may 
cause certain initial conditions to develop singularities in finite time, since this is precisely 
what happens for simpler operators. 

We consider the following simplification of T5: 

Q^ = ]^{x'D^ + D^ox^). 

One may explicitly find the deficiency subspaces, the Hamilton fiow, and the propagator 
(which may be viewed as a Fourier integral operator). And there exist smooth L^(]R) 
initial conditions, namely, 

Uq{x) = (x)"", i < a < 1, 

which develop singularities, traveling along the Hamilton fiow. 
For an example closer yet to Xs, one may consider 

Q2 = ^{x^hD, + hD, o x^) - 5h^D^ 
= x^hDx — ihx — bh^Dx- 
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For reference, we note that the propagator is 

_itQ2 / . , / \ 1 / f /'x cosh(at) + a sinh(at) \ \ 

e- . Uo{t,x) = a(xsmh(at) +acosh(at))- u, (^ ^ .^.^ (at) + a cosh (at) J J 

where a = \/bh. 

To study creation of singularities in general, one might use the methods of Jared Wunsch 
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Appendix A. The Weierstrass ^-Function 

In Section |VT] the Weierstrass p-function played a central role, as elliptic functions were 
found to parametrize the (/i-dependent) Hamilton flow of 

1 3 
Po(a;, ^; h) = -x{x^ + ^^) - -hx. 

Here we briefly describe some useful facts about the p-function, in the context of our 
problem. In what follows, we occasionally state (but do not prove) standard results, 
taken from the textbooks of Ahlfors |1] and Koblitz |14j . 

The Weierstrass p-function is a doubly periodic meromorphic function on the complex 
plane, say, with periods u}i,uj2 G C, with a double pole at each point of the period lattice, 
including the origin. And it is a solution of the differential equation 

(25) [p\z)f = A[p{z)f - g,p{z) - g,, 

where the "invariants" g2 and g^ characterize p just as well as ui and uj2. In our case, 

92 = -^^ and 03 = -Cn ■ 

Here Cq := po{x,^;h) is conserved under the flow. Since the differential equation has 
constant coefficients, p{z + zq) is also a solution, for any zq G C. 
One may factorize the cubic polynomial ( l25l) : 

p'{zf = 4{p{z) - ei){p{z) - e2){p{z) - eg). 

It is a familiar fact that the roots 61,62, 63 of 4x^ — g2X — g^ are distinct if and only if the 
discriminant is nonzero: 

A := gl - 21gl ^ 0. 

Moreover, all the roots ei, 62, 63 are real if and only if g2 and g^ are real and A > 0. In 
our case. 



28 



MICHAEL VANVALKENBURGH 



So if we ignore the trivial case Cq = (see Section rVTj) . we have that 61,62, 63 are distinct 
if and only if C^^h^, and moreover 61,62, 63 are real if and only if Cq < h^. 

Suppose that we are in the case A > 0. We then order the Cj such that 61 < 63 < 62- 
Then one may choose the periods ui and 002 of p to be given by 

in dt 
(26) -uji = ■ 



^3 + 92t - 4t3 
dt 



and 

where we take the positive branch of the square root and integrate along the real axis. 
Hence the fundamental domain of p in the complex plane is a rectangle with its opposite 
edges identified. 

In our case, in Section IVll it is clear that the unbounded flow lines are described by 
restricting p to the real line, since the poles lie on the real line. Also, it is clear that the 
nondegenerate flow loops are described by restricting p to a line of the form 

La = {t + a; t E M}, for some a G C. 

To determine the value of a, we make use of the addition formula 



p{z + u) = -p{z) - p{u) + 



1 fp'iz)~p'iu) 



2 



ei = , 62 = py—], 63 = p 



P = 0, P Hr = 0, and p ' 



p{z) - p{u) 

For p(t + a) to be real for all t G M, we then only need a to satisfy p{a) G M and 
p'{a) G M. But it is well known that 

iUi\ f^2\ f iUi + UJ2 

and that 

2 J ' ° V 2 

Hence, in light of (1261) and (1271) . we may simply take a = ^Ui. 

This shows that the nondegenerate flow loops are parametrized by 

R3t^ix{t),^{t)), 



with 



and 



^it) = \C, (p (t + y ) - \h 



-P' [t + f ) 



This can easily be seen in the figures. The fundamental domain of p is pictured in 
Figure El As t traverses the horizonal line (a), {x{t),^{t)) traces out the curve (a) in 
Figure [31 And as t + ^ traverses the horizontal line (6), {x{t),^{t)) traces out the curve 
(6) in Figure [3l The parameters in Figure [3l are = 0.1 and Co = 0.025. Note in 
particular that Cg < h? . When Cq = h?, the loop reduces to a point, although there 
is still an unbounded component. And there is only an unbounded component when 
Cl > h\ 
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Figure 2. The fundamental domain of p. When restricted to either the 
hne (a) or the hne (6), p is real- valued. 





b 


2 
1 




a 

1 




-0.6 






0.2 


0.4 








-1 












-2 









Figure 3. The elliptic curve parametrized by {x{t),^{t)). The parameters 
for the curve are h — 0.1 and Co = 0.025. 
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